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Abstract 

TORCH is a time-of-flight detector concept using Cherenkov light to provide charged particle identification 
up to lOGeV/c. The concept and design are described and performance in simulation is quantified. 
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1. Introduction 

Charged particle identification is critical for a 
heavy flavour physics experiment. It is required 
both for separating signal decays from more co- 
pious backgrounds, such as B° s — > D~K + from 
— > D~ir + , and for tagging the initial flavour of 
neutral mesons which may undergo mixing. One 
way to accomplish this is with a Ring Imaging 
Cherenkov (RICH) system, in which mass hypothe- 
ses are tested by measuring the angle (0c) and num- 
ber of Cherenkov photons emitted by a charged 
particle passing through a radiator. For a given 
detector geometry, the separation power between 
K and 7r hypotheses depends on the refractive in- 
dex n of the radiator and the momentum p of the 
particle. Between the pion and kaon thresholds, 
pions are positively identified but kaons are iden- 
tified only by the absence of photons. In principle 
this can provide good separation, but it is vulner- 
able to background from fake tracks (e.g. due to 
mistakes in pattern recognition in the tracking sys- 
tem), which cannot create Cherenkov photons and 
so would be identified as kaons. Likewise, protons 
would not be distinguished from kaons. Therefore, 
it is desirable to work above the kaon threshold — 
and in an environment with high occupancy in the 
tracking system or a large number of protons, it is 
essential. 

At LHCb particle identification is needed for 
the momentum range from 2GeV/c to lOOGeV/c. 
This is currently accomplished by a RICH sys- 
tem Q consisting of two detectors with three ra- 
diators. Two of these are gaseous: C4F10 and CF4, 



with refractive indices of 1.0014 and 1.0005, respec- 
tively, for visible light at STP. These correspond to 
kaon Cherenkov thresholds of 9.3 and 15.6GeV/c, 
and together provide better than 3a kaon-pion sep- 
aration up to lOOGeV/c. Positive kaon identifica- 
tion from 2-10 GeV/c requires a third radiator with 
a kaon Cherenkov threshold around 2GeV/c, cor- 
responding to n m 1.03. This excludes conventional 
materials: the phase transition from gas to liquid or 
solid is associated with a large jump in (n— 1) from 
0(1O~ 3 ) or smaller to (9(1). The low density of 
aerogel gives it a suitable refractive index: LHCb's 
third radiator is aerogel with n — 1.029 in air and 
n = 1.037 in C4F10. However, Rayleigh scattering 
limits the useful light yield: approximately 7 de- 
tected photoelectrons are expected at LHCb for a 
saturated track in aerogel. In a high-background 
scenario, such as the proposed LHCb upgrade from 
2 x 10 32 cm _2 s _1 to 2 x 10 33 cm~ 2 s _1 , kaon-pion 
discrimination from aerogel would be severely com- 
promised 

Another method of charged particle identification 
is by time of flight measurement. The principle 
is simple: if a particle of mass m is detected at 
position x and time t after leaving the origin, then 
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and so the difference in expected time between 
kaons and pions would be 
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Therefore, to obtain 3cr separation between the two 
hypotheses the time resolution at must satisfy 
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If this approach were to be used for kaon iden- 
tification in the range 2-10GeV/c and the path 
length were x — 10 m, then a per-track time res- 
olution (7 t < 12.5 ps would be required. A very 
fast response would be needed — Cherenkov radia- 
tion would be suitable for this purpose. If the ex- 
pected light yield were 50 detected photoelectrons, 
a per-photon time resolution under 90 ps would be 
required. 

A third method is to measure the time of propa- 
gation t 7 over a path length d 7 of Cherenkov pho- 
tons emitted by a charged particle. In a dispersive 
medium, photons propagate at the group velocity 
c/n g , and so 
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The relationship between n and n g is a non-linear 
function that depends on the medium but obeys: 
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After computing n from n g , the mass of the particle 
can be extracted: 
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The methods outlined above may be combined 
in an elegant way. Consider a large, thin plane of 
an optically dense medium such as quartz. Charged 
particles passing through the plane emit Cherenkov 
photons, which propagate to the edges via total in- 
ternal reflection as in the BABAR DIRC 0]. The 
arrival time of a photon at the edge of the plane 
is the sum of the time of flight of the track to the 
emission point, i, and the time of propagation of 
the photon in the medium, r 7 . Then 




where n g may be determined from 
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Both dt/dm and dr 7 /dm have the same sign so 
the two effects combine constructively. For a given 
track, the first term of Eq. [5] is fixed but the sec- 
ond depends on the direction and path length of the 
photon, which must be reconstructed for each pho- 
ton individually — for the planar geometry outlined 
above, this can be done by measuring two indepen- 
dent angles for each photon. This is the principle 
behind the design of the TORCH, and of the Belle 
upgrade TOP [5] and PANDA endcap DIRC [1] 
which inspired it. 

2. TORCH design 

The TORCH detector consists of a plane of 
quartz covering the full LHCb angular acceptance 
and placed 12 m downstream of the interaction 
point, plus focusing elements and photodetectors 
along each edge. The quartz plate has dimensions 
7440 mm x 6120 mm, with a hole of 26 mm x 26 mm 
at the centre for the beampipe, as illustrated in 
Fig.[l] The focusing element uses a cylindrical mir- 
ror to focus the light such that the distance along 
the photodetector plane is approximately propor- 
tional to the angle 9 Z , as illustrated in Fig.[2j This 
angle is key: the overall resolution is dominated by 
the pixelization uncertainty on 9 Z . From simulation 
studies, we find that a 128-channel segmentation 
and a dynamic range of 400 mrad corresponds to 
an uncertainty of 0.96 mrad on 9 Z and per-photon 
uncertainty on the time of propagation of O(70ps). 
The granularity requirement in the other dimension 
is much looser, since the azimuthal angle is typically 
measured with a lever arm of a few metres: 0(1 cm) 
is sufficient match the 1 mrad resolution. 

The key property required of the photodetector 
is high speed: its intrinsic time spread should be 
significantly better than 90 ps. The photon detec- 
tion technology which currently gives the best time 
resolution is the use of micro-channel plate photo- 
multipliers (MCP-PMT). They have been used in 
R&D for a related concept under development by 
the Belle collaboration for their upgrade [7], with 
single-photon time resolutions of 35 ps achieved. 

An example of a commercially-available detector 
of this type which comes close to satisfying the re- 
quirements of TORCH is the Photonis™ XP85022 
MCP-PMT. This is available in a 1024-channel ver- 
sion, with an array of 32 x 32 pixels in a 53 x 53 mm 2 
active area, and overall dimension of 59 mm. For 
the sake of developing a concrete design for TORCH 
we have assumed the use of this photodetector, laid 
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Figure 1: Schematic layout of the TORCH detector. For 
clarity, focusing elements have only been shown on the upper 
and lower edges. Photodetectors extend along the full length 
of each edge. 
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Figure 2: Cross section of the focusing element (left), and 
definition of the angle 9 Z (right). The focusing of photons is 
shown for five illustrative angles between 450 and 850 mrad, 
emerging at different points across the edge of the plate. 

side-by-side with a pitch of 60 mm along the edges 
of the plane. The anode structure of an MCP-PMT 
consists simply of conductive pads to read out the 
charge, and the pixcllization can therefore in prin- 
ciple be adjusted. For the simulation of TORCH we 
have assumed an array of 8 x 128 pixels of dimension 
6.6 mm x 0.4 mm. 

3. Reconstruction and pattern recognition 

Following the approach outlined in Sec. [1] for the 
design from Sec. [21 the reconstruction is relatively 
straightforward: we measure the arrival time and 
position of a photon on the photodetector plane, in- 
fer 9 Z , and calculate its trajectory assuming it was 
emitted by the track at the midpoint in z of its path 



through the quartz block. Given measurements 
of the track path and momentum from a tracking 
system and knowledge of the optical properties of 
quartz, this is sufficient information to compute all 
quantities in Eq.[5]and[§]and extract the mass from 
t + r~ f . There are two major practical challenges, 
though: first, association of photons to tracks (pat- 
tern recognition); and second, determining the time 
when the track left the origin. 

Pattern recognition is critical: in the environ- 
ment of the upgraded LHCb luminosity, there may 
be 0(100) fully reconstructed tracks plus a large 
number of secondaries passing through the TORCH 
in an event. However, most photon-track pairs can 
be rejected as unphysical. We take advantage of 
the limited range of photon wavelengths to which 
the photodetector is sensitive and reject candidates 
outside this range. This is mathematically equiva- 
lent to a restriction on 6c '■ for our geometry, this 
means that relevant photons lie roughly along arcs 
on the photodetector plane. Additionally, we can 
ignore photon-track pairs with unphysical timing. 
Instead of attempting to measure the mass of the 
particle directly, we assume e, [i, n, K,p hypotheses 
and test for consistency with each in turn. Back- 
ground photons whose timing is not consistent with 
any of these masses can be ignored. 

In Sec. Q] it was assumed that the track left the 
origin at t = — or equivalently, that the relative 
timing between the photodetector and the initial 
collision is known. However, we cannot rely on an 
external clock for lOps-level timing — and even if 
we could, the luminous region at LHCb extends a 
few cm in z, which would add a smearing of order 
a few tens of ps. We must therefore measure the 
track start time ourselves. One approach would 
be to install a second TORCH-like detector close 
to the interaction point, at the price of additional 
cost, material, and time resolution. A more efficient 
solution is to take advantage of the high pion mul- 
tiplicity at a hadron collider and use other tracks 
from the same event primary vertex to fix the rel- 
ative timing. The same reconstruction procedure 
described above is used except that the last step is 
inverted: for each track from the event primary ver- 
tex, the pion mass is assumed as the input and the 
time elapsed (t + T 7 ) is calculated from Eq. [51 Sub- 
tracting this from the measured photon detection 
time gives a measurement of the track start time. 
For a primary vertex with fully reconstructed 
pion tracks, the start time resolution is smaller than 
the per-track time resolution by roughly y/N^. 
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4. Performance 

To test the particle identification performance, 
we take events from the full LHCb Monte Carlo 
simulation, determine the set of charged particles 
that would pass through the TORCH, simulate the 
emission, propagation, and detection of Cherenkov 
radiation in the TORCH for these particles with a 
stand-alone program, then apply the reconstruction 
and pattern-recognition procedure outlined above 
to the output. No truth information is used for the 
reconstruction and pattern-recognition. However, a 
number of simplifying assumptions were made: per- 
fect measurements of track position and momen- 
tum, perfect extrapolation of tracks through the 
magnetic field of the detector, negligible electronics 
noise and cross-talk; negligible spill-over between 
events; and no multiple scattering, delta ray emis- 
sion, or inelastic collisions in the TORCH. Back- 
ground from certain sources such as backscatter 
from the calorimeters was not available to the sim- 
ulation and was also neglected. 

We characterize the performance in terms of 
particle identification efficiencies for those kaon 
and pion tracks which would be useful for physics 
analysis. Specifically, we require tracks to be 
well-measured (with track segments found both 
upstream and downstream of the magnet), well- 
matched to an event primary vertex, and associated 
to a kaon or pion in the Monte Carlo truth infor- 
mation. Tracks are identified either as kaons or as 
pions according to which of the two hypotheses is 
found to be most consistent by the pattern recogni- 
tion. The efficiency is shown as a function of track 
momentum in Fig. [3] for events simulated at a lu- 
minosity of L = 2 x 10 32 cm _2 s _1 . The efficiency 
for correct identification is > 95% up to lOGeV/c, 
dropping at higher momentum as the time differ- 
ence between kaon and pion hypotheses diminishes. 
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Figure 3: Identification efficiency of the subset of well- 
measured charged tracks which are well-matched to a pri- 
mary vertex. The plots show the efficiency for a kaon (up- 
per) or pion (lower) track to be identified correctly (black) 
or incorrectly (red/dot), considering only the kaon and pion 
hypotheses. 



modular design with a smaller quartz plane is also 
under investigation. 
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5. Conclusions 

The TORCH concept is intended to provide 
charged particle identification in the momentum 
range 2-10GeV/c in a high-rate environment. It 
has been shown that the pattern-recognition prob- 
lem is solvable and that target performance can 
be reached under simplified conditions. The chal- 
lenge is now to test whether this performance can 
be maintained under realistic conditions and at 
L = 2 x 10 33 cm~ 2 s _1 . This will require proto- 
type tests and more detailed simulation. A revised, 
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